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Q-r We present a weak lensing analysis of the X-ray luminous cluster CL 0152- 

■ 1357 at z ~ 0.84 using HST/ACS observations. The unparalleled resolution 

■ and sensitivity of ACS enable us to measure weakly distorted, faint background 
galaxies to the extent that the number density reaches ~ 175 arcmin -2 . The PSF 
of ACS has a complicated shape that also varies across the field. We construct a 

^ ■ PSF model for ACS from an extensive investigation of 47 Tuc stars in a modestly 

crowded region. We show that this model PSF excellently describes the PSF 
q\ ■ variation pattern in the cluster observation when a slight adjustment of ellipticity 

is applied. The high number density of source galaxies and the accurate removal 
of the PSF effect through moment-based deconvolution allow us to restore the 
^ . dark matter distribution of the cluster in great detail. 
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The direct comparison of the mass map with the X-ray morphology from 
Chandra observations shows that the two peaks of intracluster medium traced 
by X-ray emission are lagging behind the corresponding dark matter clumps, 
indicative of an on-going merger. The overall mass profile of the cluster can 



be well described by an NFW profile with a scale radius of r s = 309 ± 45 kpc 
and a concentration parameter of c = 3.7 ± 0.5. The mass estimates from the 
lensing analysis are consistent with those from X-ray and Sunyaev-Zeldovich 
analyses. The predicted velocity dispersion is also in good agreement with the 
spectroscopic measurement from VLT observations. In the adopted cosmology 
where f2 M = 0.27, f2 A = 0.73, and h = 0.71, the total projected mass and the 
mass-to-light ratio within 1 Mpc are estimated to be (4.92 ± 0.44)xl0 14 M Q and 
95 ± 8 M q /Lbq, respectively. 
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1357) — galaxies: high-redshift 

1. INTRODUCTION 

Gravitational lensing has been a unique tool to probe the intervening matter distribution 
between the observer and source objects without any assumption about the dynamical phase 
of deflectors. The most impressive and beautiful images of giant arcs can appear when the 
light from background sources passes near the caustic of a foreground lens. These "strong 
lensing" features indicate the presence of a critical (or higher) surface mass density in the 
inner region and are particularly useful for constraining the core structure of the lensing 
cluster. In the "weak lensing" regime where the distortion becomes weaker and less obvious 
(whether the cluster is less massive or projected distances of source galaxies are farther 
from the cluster center), one can still detect coherent alignments of background galaxies 
and restore the mass distribution up to far greater radii from these subtle measurements. 
Since the first successful detection of systematic alignments of background galaxies by Tyson, 
Wenk, & Valdes (1990), the technique has been applied to a wide selection of galaxy clusters 
and is now firmly established as one of the most straightforward paths to probe the mass 
distribution of the cluster in question. 

In general, the detectability of weak gravitational shears depends not only upon the 
intrinsic signal strength determined by the projected mass density of a lens and the geometry 
between lens and source, but also upon the observational restrictions set by finite sensitivity 
and resolution of an instrument. In this regard, a weak lensing analysis of high-redshift 
clusters is disadvantaged because the signal decreases as the redshift of the lensing cluster 
approaches that of background galaxies, and also it becomes harder to recover shapes of faint, 
poorly resolved galaxies, which however contain most of the useful signal. Nevertheless, the 
demands for comprehensive studies on many individual high-redshift clusters are increasing 
because of their potentially significant implications for cluster formation and cosmology. For 
example, even the mere abundance of massive clusters at such high redshifts can strongly 
constrain the cosmological density parameter Q, decoupling the agQ ' 5 — 0.5 degeneracy 
(e.g., Carlberg, Morris, Yee, & Ellingson 1997; Bahcall & Fan 1998). Furthermore, most 
high-redshift clusters possess filamentary structures indicating their early stage of formation, 
and it is interesting to investigate in detail how dark matter is distributed with respect to 
cluster galaxies or the intracluster medium (ICM) traced by X-ray emission. 

The remarkable substructure of MS 1054-03 obtained through the weak lensing analysis 
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of WFPC2 observations (Hoekstra, Franx, & Kuijken 2000, hereafter HFK00) demonstrated 
the undeniable merits of space-based weak lensing observations and already hinted at the 
bright prospects of the newly installed Advanced Camera for Surveys (ACS) in this ap- 
plication. The advantages of HST observations over ground-based imaging include a higher 
number density of background galaxies whose shapes can be reliably determined, and smaller 
corrections for point spread function (PSF) effects. The pixel size of the Wide Field Channel 
of ACS is ~0.05 ", offering twice the sampling resolution of the Wide Field (WF) chips of 
WFPC2. In addition, ACS has a factor of 5 improvement in throughput while providing 
twice the field of view of WFPC2. Therefore, the finer resolution, higher sensitivity, and 
wider field of view of ACS can provide a much higher fraction of well-resolved galaxies with 
less investment of HST observing time. 

In this paper, we present a weak lensing study of the X-ray selected high-redshift 
(z-0.84) cluster CL 0152-1357 using the Wide Field Channel of ACS. Together with MS 
1054-03, CL 0152-1357 is one of the most X-ray luminous clusters at z > 0.8 known to 
date whose X-ray properties have been well-studied by the ROSAT (Ebeling et al. 2000), 
BeppoSAX (Delia Ceca, Maccacaro, Rosati, & Braito 2000), and Chandra Observatory 
(Maughan et al. 2003). Nevertheless, unlike MS 1054-03, there have been no HST-based 
high-resolution weak lensing studies so far. Such studies, combined with X-ray observations, 
can substantially enhance our understanding of the dynamical evolution of the ICM and its 
interaction with cluster galaxies as well as of the cluster substructure as a whole. Particu- 
larly, CL 0152-1357 is known to have a complicated X-ray and optical substructure, and the 
high-resolution mass reconstruction from the ACS weak lensing is expected to provide the 
mass distribution of the cluster in unprecedented detail. 

We organize our works as follows. In §2, we describe the observations and basic data 
reduction. Ellipticity measurement of source galaxies is presented in §3. In §4 we discuss the 
PSF modeling and correction. §5 describes the luminosity of the cluster and the distribution 
thereof. The mass reconstruction of the cluster is handled in §6 and the total mass estimates 
from various approaches are described in §7. Finally, the substructure from this study is 
compared with that of other studies in §8 before the conclusion §9. Throughout the work we 
assume A CDM cosmology favored by the Wilkinson Microwave Anisotropy Probe (WMAP) 
with VL M = 0.27, fi A = 0.73, and H = 71km s _1 Mpc _1 . 



2. OBSERVATIONS 



CL 0152-1357 was observed in a 2x2 mosaic pattern allowing ~ 50" overlap between 
pointings with the Wide Field Channel of ACS during 2002 November and 2002 December 
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(GTO proposal 9290, P.I. Ford). The cluster was imaged in F625W, F775W, and F850LP 
(hereafter ^25,^775, and z 850 , respectively) with integrated exposure per pointing ~ 4800 s. 
The low level CCD processing (e.g. overscan, bias, dark subtraction, and flat-fielding) was 
done with the standard STScI CALACS pipeline, and the final high-level science images 
were created through the "apsis" ACS GTO pipeline (Blakeslee et al. 2003). The result- 
ing high accuracy in both the image registration and the geometric distortion correction 
is indispensable in weak lensing measurements. The apsis pipeline calculates offsets and 
rotations between images using the "match" program (Richmond 2002) after applying a 
geometric distortion model (Meurer et al. 2003) to the astronomical objects which have a 
high signal-to-noise ratio (SNR). The typical shift uncertainty is ~ 0.015 pixels (J. Blakeslee 
2004, private communication). The drizzle-blot-drizzle cycle of apsis automatically rejects 
cosmic rays and generates mosaic science images. We used the Lanczos3 drizzling kernel, 
which gives a sharper PSF and less noise correlation between neighboring pixels. The noise 
correlation decreases the root mean square (RMS) noise fluctuations and therefore causes 
photometric errors to be slightly underestimated. Apsis calculates the correct RMS noise in 
the absence of correlation. The RMS map produced in this way is used for source detec- 
tion and photometric error estimation. We present the color composite image of the cluster 
center in Figure 1, which shows many strong-lensing features such as arc(let)s around the 
two central brightest cluster galaxies (BCGs). The objects were detected using SExtractor 
(Bertin & Arnouts 1996) on the detection image created by combining all the present filter 
images after applying inverse variance weighting. In general, many parameters in SExtractor 
affect the detection procedure, and we adopted the values obtained from the experiments of 
Benitez et al. (2004). They tried to suppress spurious detections while extracting all obvious 
galaxies by tuning up these parameters. We found that their rather conservative choice of 
DETECT_M I NAREA= 5 and DETECT_THRESH= 1 . 5 selects faint galaxies up to the detection limit 
(S/N ~ 3) without significantly introducing false objects into the object catalog. We ran 
SExtractor in dual image mode for each filter to obtain the galaxy photometry and colors. 
The final catalogs were visually compared with the image in order to remove false detections 
(e.g. diffraction spikes and spurious spots around bright stars, saturated CCD bleeding, 
noise fluctuation at image edges, clipped and merged galaxies, HII regions inside nearby 
galaxies, etc.). This resulted in a total of 10992 objects. The number count plot in Figure 2 
shows the completeness of our data down to ~27.5 magnitude in all three passbands. 



3. ELLIPTICITY MEASUREMENTS 



In the regime where the background galaxy is much smaller than the scale length of 
the gravitational potential variation, we can obtain the linearized lens mapping equation as 
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follows: 



A(x) = 5., 
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where A(cc) is the transformation matrix x' = Aa; which relates a position x in source plane 
to a position a;' in image plane, and \& is the 2-dimensional lensing potential. In the matrix 
of equation 1, the convergence k determines the overall magnification, and 71 and 72 describe 



the shear along x-axis and at 45 from the x-axis, respectively. Therefore, in general, the 
galaxy images are distorted in shape (71, 72) and size («) (in a strict sense, 71(2) also alters 
the object size anisotropically) . Though Broadhurst, Taylor, & Peacock (1995) claimed that 
the magnification bias (and thus the number density bias) can be used to estimate the local 
surface mass density directly, most weak lensing works have been based on the ellipticity 
biases. This is because the magnification effect is more sensitive to shot noise, and the SNR 
of the shear measurement is considerably better than that of the magnification effect in a 
typical weak lensing (k <C 1) regime (Schneider, King, & Erben 2000). 

The ellipticity of an object can be defined in terms of weighted quadrupole moments as 
follows: 



where f(x) is the pixel intensity, and w(x) is the weight function required to suppress the 
noise in the outer region of the object. Kaiser, Squires, & Broadhurst (1995, hereafter KSB) 
used the circular Gaussian weight whose size matches that of the object, thus maximizing the 
significance of the measurement. However, the circular weight makes the object rounder, and 
the effect becomes severe for highly non-circular galaxies. Besides, the ellipticities calculated 
in this way do not follow the simple ellipticity transformation rule (Kochanek 1990; Miralda- 
Escude 1991) in response to the applied shear. These features necessitate the introduction 
of an additional parameter which in general depends on the higher moments. In KSB work, 
this quantity is referred to as "shear polarizability" P 1 . 

Recently, Bernstein & Jarvis (2002, hereafter BJ02) introduced adaptive moments using 
an elliptical weight function whose shape and size match those of an object. While the con- 
cept of finding the optimal elliptical Gaussian weight function is mathematically simple, the 
actual implementation can take various forms. For example, one can determine the weight 
function by minimizing the deviation from the image in the least-square sense. Alterna- 
tively, one can start with a circular weight function and iteratively modify the ellipticity, the 
size, and the centroid of the weight function until these parameters converge. BJ02 effected 
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the determination of optimal elliptical weight function by iteratively shearing the objects to 
match the circular Gaussian weight. Considering the finite pixelization of object images, 
this may not sound more attractive than the previous two schemes. However, if the galaxy 
images can be decomposed via mathematically well-behaved basis functions, the adaptive 
elliptical moments are computed inexpensively. BJ02 proposed the polar eigenfunctions of 
2-dimensional quantum harmonic oscillator (QHO) as basis functions. This decomposition 
was also independently suggested by Refregier (2003, hereafter R03) though his shear esti- 
mator is different from that of BJ02. Many mathematically convenient formalisms developed 
for these eigenfunctions or shapelets include operators which can effect coordinate transfor- 
mations such as shear, translation, dilation, rotation, etc. Even more important advantages 
one obtains from the galaxy expansion using shapelets are that the PSF can be compactly 
described by the coefficients of the basis functions and the (de) convolution is easily achieved 
by simple matrix manipulations. 

Shapelets in polar coordinate are given by 

I(r,9)= 5> p ^(r,0) (4) 

P,q>0 

where L™(x) are the Laguerre polynomials. The complex conjugate relation ty° p = ^/^ p is 
used to compute the eigenfunctions when p < q. Other useful mathematical properties of the 
above basis functions along with recursion relations of matrix elements of the aforementioned 
operators are presented in BJ02. 

Once the galaxy image is decomposed into b pq vectors, we can translate, dilate, and 
shear 

b' = (S V V^T Z ) ■ b, (6) 
until the series of transformation satisfies the following conditions 
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The condition imposed by equations 7, 8, and 9 relate to centroid, size, and ellipticity of 
the optimal elliptical Gaussian, respectively. In this paper, the algorithm of BJ02 method is 
independently implemented in the Interactive Data Language (IDL). Figure 3 demonstrates 
the statistics of the decomposition of galaxies in r 62 5 passband. The fraction of reliably 
measurable galaxies decreases substantially if target galaxies are fainter than r$25 ~ 29. A 
similar trend is observed in the other two passbands. The final shape catalog is produced 
after optimally combining ellipticities in all three passbands. 

One must note that the quantity n in equation 6 is different from the conventional 
definitions of ellipticities though they are related in a straightforward manner. 

1 - a 2 

6 = Y+lf =tanhr ? ( 10 ) 

e =I^i = tanh^, (11) 

where q is the axis ratio b/a. In this paper, the distortion 5 in equation 10 is referred to as 
ellipticity unless indicated otherwise. 



4. PSF CORRECTIONS 

As one probes to weaker and weaker lensing regimes, the accurate removal of any instru- 
mental artifact becomes paramount to the success of the analysis. Finite seeing causes the 
cicularization of small galaxies while the anisotropy of the PSF can create systematic biases 
in the size and the direction of the polarization. The pioneering investigation of KSB95 sug- 
gested that an approximation can be made by treating the real PSF as a small perturbation 
to an isotropic PSF. The original prescription and its variations have been employed widely 
during the last decade though the valid regime of its application was frequently questioned. 
Kaiser (2000) argued that despite the fact that the modified KSB works reasonably in some 
cases, the technique may become problematic in diffraction-limited observations. In the cur- 
rent investigation, among many recent suggestions (e.g., Wilson, Cole, & Frenk 1996; Kuijken 
1999; Kaiser 2000), we settled upon the moment-based deconvolution technique (BJ02; R03) 
which performs the deconvolution by the matrix manipulation of shapelet components of 
the galaxy and PSF. The decomposition of the PSF in this way not only eases the modeling 
of the PSF variation across the field, but also effectively suppresses the noise amplication if 
the truncation of the higher order moments is carefully handled. 
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4.1. PSF Modeling of WFC 



Though field-dependent variation of the WFC PSF is small compared to that of WFPC2, 
its change in ellipticity within the field is significant (Krist 2003). In order to investigate the 
issue we retrieved the repeated (~ every three weeks) observations of the modestly crowded 
region of the globular cluster 47 Tuc originally used to monitor the flat-fielding stability 
of ACS (PROP 9656, PI De Marchi). Because the default drizzling kernel of the STScI 
pipeline is square, we had to re-drizzle all the flat-fielded (FLT) images using Lanzcos3 
kernel to match the PSF size of the CL 0152-1357 observation. After the initial detection 
of stars by SExtractor, we selected "good" stars which are bright, unsaturated, and isolated 
(having no companion stars or cosmic-rays within 15 pixels from the center). Then, each 
star is decomposed into shapelet components by minimizing the following: 



where the optimal centroid and size of the eigenfunctions ^> pq are iteratively determined. 
The two whisker plots in Figure 4 show typical PSF patterns of the WFC measured from 
the 47 Tuc field. Krist (2003) pointed out the magnitude of ellipticity is determined by the 
focus offset, and the angle of elongation switches by 90° when the sign of the offset becomes 
opposite. This fact is confirmed by comparing Figure 4a and 4b which are taken on 2002 
October 3 and 2002 October 24, respectively, roughly on the same patch of the 47 Tuc field. 
We find that ACS PSF patterns at other epochs follow one of these two patterns with slightly 
altered ellipticities. 

The spatial variation of shapelet coefficients b pq of ACS PSFs is modeled as 



We found the third order in x t y j (i.e. i+j < 3) is sufficient to describe the pattern and 
higher order polynomials do not improve the agreement between the model and the data. 
Now, the important question is how well the PSF taken from the 47 Tuc field can describe 
the PSF on the CL 0152-1357 field. Since the ellipticity of the PSF changes with respect to 
the focus offset, it is necessary that the ellipticity of the model PSF is made adjustable to 
match that of the actual PSF in the cluster observation. We implemented this by applying 
a shear operator to the shapelet components of the model PSF. That is, 
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where the evaluation of matrix elements of the shear operator Ss v can be found in Appendix 
A. 3 of BJ02. Though Sr] can be also allowed to vary depending on the position in principle, 
we found that a simple fixed parameter per exposure nicely reduce the systematic residuals. 
Because we measure galaxy shapes on the 2x2 mosaic image, another slight complexity 
arises due to the overlap between pointings. Nevertheless, a reasonable assumption can be 
made that the PSF in the overlapping region is very closely approximated as an exposure- 
time-weighted average of all the contributing PSFs. 

Using a typical half-light radius versus magnitude plot, we initially selected 73 isolated, 
bright stars which can be used as local PSF indicators. We removed stars having any 
noticeable defects from the list by visual inspection, which ended up a total of 62 stars. 
Figure 5a shows the polarization pattern in CL 0152-1357 observation measured from these 
stars. Then, using our model PSF we constructed "rounding kernels" (Fischer & Tyson 1997; 
Kaiser 2000; Bernstein & Jarvis 2002) which circularize the originally elongated PSFs and 
applied them to the CL 0152-1357 images. Comparing the ellipticities before (Figure 5a) and 
after (Figure 5b) the application of the rounding kernel verifies that our model PSF closely 
represents the real PSFs on the cluster image (see also Figure 6). It appears that there still 
remain tiny but systematic residuals due to the incompleteness of the model; however, their 
effects on the cluster mass analysis are estimated to be negligible. 



Though the "rounded images" obtained in the previous section can be used to mea- 
sure the object shapes, this is not preferred to the straightforward deconvolution technique 
because of the following reasons. First, the rounding kernel always degrades the original 
image seeing because the kernel size must be comparable to the instrument PSF size in 
order to remove the anisotropy sufficiently, which in particular is detrimental to very small 
galaxy images. Second, the dilution (circularization) correction provided by BJ02 is still 
an approximation and Hirata & Seljak (2003) showed in their simulation that indeed the 
prescription by BJ02 is not accurate if the kurtosis of the PSF is not small or the galaxy is 
not well-resolved. 

Due to the simple transformation rule of the Gaussian functions, the deconvolution can 
be effected by convenient matrix manipulations: 



4.2. 



PSF Correction from Deconvolution 




(15) 
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where b°, b\ and b s are shapelet components of the convolved image, the pre-seeing image, 
and the PSF, respectively. 

The evaluation of the matrix elements Cl iq n iPaqa is summarized in BJ02 (see also R03 for 
Cartesian coordinates). After contracting C , £** Ps9s and b s Psqs , we get 

b° = V pw«6< . (16) 
Now P Pi ^ can be inverted to compute the deconvolved image &!, „. from the PSF-convolved 

Polo 1 ^ PiHi 

original image b° „ . Because we desire to make the matrix P Pi f invertible and also minimize 

Po<4o Po^lo 

the noise amplification which is typical in every deconvolution problem, the expansion of the 
PSF in shapelets must be truncated appropriately and the characteristic size of the object 
should be large enough compared to the size of the PSF. 



5. LUMINOSITY ESTIMATION 

We base our selection of cluster members on the tight color-magnitude (CM) relation 
of early-type galaxies of the cluster. Because the 4000A break of the cluster ellipticals 
is redshifted to the cutoff wavelength of r$25 filter, (tq25 — -285o) colors are better suited 
than («775 — -2g5o) colors. As shown in Figure 7, the bright cluster red sequence of CL 
0152-1357 occupies a relatively narrow strip in the (r 62 5 — -Zsso) versus z 850 CM diagram. 
Because increasing photometric errors at faint magnitudes cause the distinction to become 
less apparent, we selected 371 galaxies brighter than z 850 ~ 25. The spectroscopic catalog 
from VLT observations (R. Demarco et al. 2004, in preparation) is used to reject bright 
non-cluster members (z 850 < 22) and to include some known blue cluster galaxies. We 
show the smoothed cluster light distribution from these member galaxies in Figure 8. The 
spectroscopic survey of the CL 0152-1357 field serendipitously discovered a foreground (z ~ 
0.63) group of ~ 12 galaxies rather loosely scattered over the entire field. We excluded these 
galaxies in the above light distribution. The vertically elongated main structure as well as 
the less luminous but distinct clumps around the main body is clearly visible. We refer to 
the brightest concentration in the light distribution as the cluster center hereafter. This 
smoothed light distribution will be compared with those of the Chandra X-ray and the weak 
lensing mass in §8. 

In order to estimate the rest-frame luminosity of the cluster, we proceed as follows. 
From the dust maps of Schlegel, Finkbeiner, & Davis (1998), we obtained E(B-V)=0.014 and 
determined the extinction corrections for i 775 and z 850 to be 0.028 and 0.020, respectively. 
Then, synthetic photometry is performed by combining the latest ACS throughput curves 
(M. Sirianni et al. 2004, in preparation) and the Kinney-Calzetti spectral templates (Kinney 
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et al. 1996) so as to establish the photometric transformation of i 775 at z ~ 0.84 to the rest 
frame B magnitude. The linear best-fit result has the following form. 

Brest = *775 - (0.39 ± 0.03) (2775 - ^ 850 ) + (0.76 ± 0.05) - DM, (17) 

where DM is the distance modulus of 43.63 in this cosmology and the uncertainties of the 
coefficients are estimated assuming an accuracy of ~2% in the synthetic photometry. The 
total luminosty of the cluster is estimated by 

L B = J2 10 OA{MB °- Brest - j) L Be , (18) 

3 

where Mb® = 5.48 is the absolute B magnitude of the sun. 

However, the L B obtained in this way does not include the contribution from the faint 
(^850 > 25) population. In addition, we expect the blue cluster members also comprise a 
significant fraction of the total light because CL 0152-1357 is a high-redshift cluster. We 
choose to adopt the scheme by HFK00 in order to correct the total luminosity of the cluster 
for these factors. By fitting the Schechter luminosity function to our sample galaxies, we 
found ~ 4% of the total light must be added in order to account for the faint population. 
To estimate the fraction of the blue cluster galaxies, we compared the spectroscopic catalog 
with the CM selection and determined that we would lose ~ 20 % of the total luminosity if 
not including the blue population. The fraction estimated here for the cluster CL 0152-1357 
is slightly higher than the value for MS 1054-03 obtained by HFK00 who quoted 16%. We 
present the cumulative light profile as a function of the radius from the cluster center in 
Figure 9. We observe that the profile becomes marginally steeper as the radius approaches 
the field edges because of the increasing contribution from blue cluster galaxies. The light 
profile, when reproduced without including the spectroscopically confirmed blue population, 
showed no such trend. In the WMAP cosmology, 1 Mpc corresponds to ~ 131" at z — 0.84 
and the total B-band luminosity within this aperture becomes 5.2 x 10 12 L BQ . 



6. MASS RECONSTRUCTION 

6.1. Source Galaxy Selection 

Careful selection of background galaxies must be made in both colors and magnitudes 
in order to maximize the available signals. We rely on the tight CM relation of the cluster to 
separate the early-type cluster members from background galaxies. Galaxies whose (r 62 5 — 
^85o) colors are bluer than (— 0.115z 8 5o + 4.2) are chosen. The distinction is not apparent 
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at -2850 ^ 26. In an idealized observation where most of galaxies are well resolved, one can 
desire to include as faint galaxies as possible because the distortion is greater for higher 
redshift objects. However, these faint galaxies are in reality not only more susceptible to 
measurement noises, but also the amount of correction due to the circularization of the PSF 
is greater, which increases the uncertainty of ellipticities by the same factor. To establish 
the limiting magnitude for background galaxy selection, we carried out the following "shear 
recovery" test. After galaxy shapes are measured in the original field of CL 0152-1357, we 
artificially sheared the entire image by 5 % in real space. Then, the ellipticities of galaxies 
are determined once more on the distorted image, and we checked how well the applied 
shear is recovered as a function of magnitude (Figure 10). We observe that in spite of 
growing uncertainty as magnitude increases, the shear is recovered down to ~ 29.5 mag. 
Another useful experiment is to examine the strength of the tangential shear while varying 
the magnitude limit of the sample. Tangential shear is defined as 

7r = ~7i cos 20 — 72 sin 20, (19) 

where is the position angle of the object with respect to the cluster center. If no shear is 
present, the average of the tangential shear measured in the annulus around the center must 
approach or oscillate near zero. However, if the shear is strong enough to be measurable, (72-) 
tends to be positive and the amplitude is proportional to the magnitude of the shear. We 
divided all the detected galaxies into "bright" (24 < z 850 < 26), "faint" (24 < z 850 < 28.5), 
and "faintest" (24 < z S50 < 30) samples, and the amplitude of azimuthal averages of the 
tangential shear for different samples are compared (Figure 11). Though it is not apparent 
whether or not the signal from "faintest" galaxies is strongest, the amplitude of tangential 
shear from "faint" and "faintest" galaxies are undeniably greater than "bright" galaxies. 

We also examined the dependence of the lensing signal on source galaxy colors by 
subdividing these "faint" and "faintest" samples into "blue" and "red" subsamples using 
r 625 — -^850 colors. We do not detect any significant change in shear strength between different 
color groups, in contrast to the result of HFK00 who reported their "blue" galaxies show 
stronger signals. However, the difference must be interpreted with the different depth of the 
observations in mind. One plausible scenario is that the faint blue galaxies (FBGs) dominate 
the relatively low redshift, brighter background population whereas the contribution from 
faint red galaxies becomes increasingly important for fainter background sources at high 
redshift s. 

Considering the results of these experiments with the stability of deconvolution, we 
choose "faint" galaxies as our "best" sample. The analysis hereafter is based solely on 
these galaxies. The average number density of source galaxies in this sample reaches ~ 
175arcmin~ 2 . 
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One of the useful methods to examine the significance of the lensing detection as well 
as the systematics is to perform the following null (cross shear) test. If the lensing signal in 
Figure 11 is arising from the gravitational lensing, the resulting shear must disappear when 
source galaxies are rotated by 45° as shown in the top panel of Figure 12. We verify that the 
amplitude of the scatters are consistent with that of the randomization test where source 
ellipticities are shuffled while galaxy positions are held fixed. The bottom panel of Figure 12 
shows the result from one realization of this radomization. 



6.2. Redshift Distribution of Source Galaxies 

Sufficient knowledge of the redshift distribution of source galaxies is essential in order 
to achieve a proper scale in subsequent discussion of mass estimates. The critical surface 
density of the cluster is given by 

c 2 

S c = (20) 

/3=(max(0,^, (21) 

where D s , Di, and D[ s are the angular diameter distance from the observer to the source, 
from the observer to the lens and from the lens to the source, respectively. Compared 
to low-redshift clusters, the critical surface mass density of high-redshift cluster is a rela- 
tively sensitive function of source redshifts. That is, the angular diameter ratio Di s /D s in 
equation 21 for high-redshift clusters changes more steeply than when the lens is at lower 
redshifts. 

It is certain that the selected source population in §6.1 actually contains blue cluster 
members as well as foreground galaxies. The presence of these non-background galaxies di- 
lute the shear signal. We attempt to estimate the fraction of the cluster galaxy contamination 
in this sample using the publicly available deep ACS images from the Great Observatories 
Origins Deep Survey (GOODS; Giavalisco et al. 2004) and the Ultra Deep Field (UDF; 
Beckwith, Somerville, & Stiavelli 2003). The comparison of magnitude distribution of our 
source sample with those from these two surveys enables us to determine the fraction of 
cluster galaxies per magnitude bin up to the cosmic variance. Because F606W (hereafter 
^606) passband is used in both surveys instead of r 62 5, the (v 60e — z 850 ) color is transformed 
to the (r 62 5 — z 85Q ) color to maintain the consistent selection criterion. At faint magnitudes 
(-^850 <) 26), the GOODS catalog is incomplete. To infer the fraction of the cluster galaxies 
in this magnitude range, we add noise to the UDF images to match the SNR of our clus- 
ter observations and detect source objects on these degraded images. It appears that the 
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contamination from the cluster galaxies is not significant over the whole magnitude range 
(Figure 13a). 

Because typical magnitudes of background galaxies are still far beyond spectroscopic 
reach, we have to estimate the mean redshift of background galaxies via photometric red- 
shift techniques. We choose to use the photometric redshift catalog of the UDF (D. Coe et 
al. 2004, in preparation) to establish the mean redshift of background galaxies in the clus- 
ter field. Combining the aforementioned ACS UDF with the NICMOS F110W and F160W 
(hereafter jno and hieo, respectively) observations, we obtain reliable photometric redshifts 
of faint galaxies well beyond the faint end [z ~ 28.5) of our source galaxy sample. The 
detailed description of the observation and the photometry 1 including the PSF matching 
across the different passbands will be presented with the public release of the photomet- 
ric redshift catalog. We generate the photometric redshift catalog of the UDF using the 
Bayesian Photometric Redshift code (Benitez 2000, hereafter BPZ). The obvious advantage 
of the BPZ over the maximum-likelihood approach includes the use of the additional infor- 
mation on probability distribution for given magnitudes termed priors. We take the redshift 
distribution of Hubble Deep Field North (HDF-N) as priors, and the spectral template li- 
braries of E, Sbc, Scd, and Im by Coleman, Wu, & Weedman (1980) are selected. We also 
added two starburst templates (SB2 and SB3) by Kinney et al. (1996). The recent post- 
launch throughput curves of ACS (M. Sirianni et al. 2004, in preparation) are incorporated 
into the synthetic photometry of these SEDs. The estimated (f3) as a function of z 850 is 
shown in Figure 13b. In order to compute the final for the cluster, however, we need to 
account for the relative number ratio of the background galaxies per magnitude bin as well 
as the cluster galaxy contamination derived above. We assume that the foreground fraction 
of the sample is similar to what we obtain from the UDF. Using the cosmological parameters 
considered in the current paper, we find {(3) = 0.282. This value corresponds to a single 
source plane at (z) ~ 1.30. 

The hypothesis that source galaxies are located in a single plane, though convenient, 
causes biases in the measurement of the reduced shear (7 = 7/(1 — k) because actual galaxies 
have a broad redshift distribution. The first order correction is approximated by (Seitz & 
Schneider 1997; HFK00) 



find that the photometric zeropoints of the NICMOS, j£f = 23.4034 and = 23.2146, released 
with the version 1.0 (2004 March 9) images still need to be refined. We adjust the NICMOS zeropoints 
based on the spectral energy distribution (SED) fitting of the galaxies whose redshifts are spectroscopically 
confirmed. 
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where g' is the uncorrected reduced shear from the single-source-plane assumption. From the 
photometric redshift catalog and equation 22, we determine that the reduced shear would 
be overestimated by (1 + 0.58k). We take into account this effect in the nonlinear mass 
reconstruction. 



Due to the intrinsic shapes of individual galaxies, the local shear at the given location 
must be estimated statistically from a population of ellipticities. Naively taking a simple 
arithmetic mean not only ignores the fact that the change in ellipticity be due to an applied 
shear 7 nonlinearly depends on the ellipticity of an object, but also discount the intrinsic 
ellipticity distribution of source galaxies. HFK00 considered the weighting scheme combin- 
ing measurement noise, intrinsic ellipticity, and preseeing shear polarizability in their weak 
lensing analysis of MS 1054-03. Kaiser (2000) and BJ02 presented a similar, but more gen- 
eralized derivation incorporating the distribution of source ellipticities. In this paper we 
adopted the "easy" weighting scheme suggested in §5.2.3 of the BJ02 paper. After smooth- 
ing the shapes of galaxies with the weighting function, we obtained the distortion field in 
Figure 14. The whisker plot obviously shows tangential alignments of background galaxies 
around the cluster center. Since the distortion computed in this way is simply related to the 
shear by 8 = 27/(1 +7 2 ) ~ 27, the shear map can be inverted to construct the parameter-free 
surface mass density of the cluster. 



In the weak lensing regime (k <C 1), the shear 7 and the dimensionless mass density k 
are simply related by the following convolution: 



6.3. 



Shear Estimation 



6.4. Weak Lensing Mass Map 




(23) 



where D(x) = — 1/(37 — ix 2 ) 2 is the convolution kernel. By simple application of the 
convolution theorem, it is straightforward to invert equation 23 to express k in terms of 7 
(Kaiser & Squires 1993, hereafter KS93), 
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K (x) = i J D*(x - x')^(x')d 2 x. (24) 

Therefore, in principle once 7 is reliably measured across the field, we can construct 
a parameter-free surface mass density map either directly from equation 24 or indirectly 
through maximum-likelihood approach where equation 23 or its variation is used as part of 
the likelihood function. The mass reconstruction obtained in this way is a useful tool to probe 
the substructure of the cluster. However, normally the direct computation of the physical 
mass from this map suffers from various artifacts of the reconstruction algorithms. Most 
serious among them are spurious negative troughs around central mass clumps, thus causing 
the total mass of the field to approach zero. This happens in part because typical variants of 
the original KS93 mass inversion algorithms perform the convolution on a finite field. But, if 
the cluster is assumed to be isolated, the boundary effect can be reduced greatly by simply 
extending the field and extrapolating the measured shear outside the original field (Schneider 
& Seitz 1995). The recent improvements such as maximum likelihood (probability) method 
(e.g., Bartelmann, Narayan, Seitz, & Schneider 1996; Squires & Kaiser 1996; Seitz, Schneider, 
& Bartelmann 1998; Marshall, Hobson, Gull, & Bridle 2002) or unbiased finite-field inversion 
(e.g., Seitz & Schneider 2001; Lombardi & Bertin 1999) can also minimize those artifacts 
without field extension. We experimented with three different algorithms: direct integration 
after smoothing (Seitz & Schneider 1995), maximum likelihood (Bartelmann, Narayan, Seitz, 
& Schneider 1996), and direct reconstruction based on variational principle (Lombardi & 
Bertin 1999). We independently implemented the first two algorithms and the last one is 
kindly provided by the authors. We refer readers to the individual paper for details of each 
method. We observe that the results from the above three algorithms show no noticeable 
differences except for minor changes near the field boundaries where the signal is expected 
to be low and biased. We present the mass reconstruction from the maximum-likelihood 
algorithm in Figure 15. 

The resemblance of this mass reconstruction to the luminosity distribution (Figure 8) 
is rather remarkable. The vertically elongated central structure is clearly visible in the 
reconstructed mass map and the locations of dominant mass peaks inside the main body 
agree well with those of cluster galaxy concentrations. We also note that even outside the 
main body the spatial correlation between the mass overdensity and luminosity peaks is 
high. We will discuss the detailed analysis and interpretation of the substructure in §8. 

Despite the high-resolution of our mass reconstruction, the surface mass map in Fig- 
ure 15 is not yet ready to be used for inferring the physical mass unless the following ambigu- 
ities are resolved. First, what we measure is the reduced shear g = 7/(1 — «) rather than the 
true shear 7. That is, we always overestimate 7 because of non-zero k. The deviation cannot 
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be ignored near the overdense regions where the assumption k <C 1 breaks down. Besides, 
the correction due to the broad redshift distribution discussed in §6.2 becomes important if 
k is not sufficiently small. Nevertheless, it is possible to correct these effects by iteratively 
updating k and 7 (Schneider & Seitz 1995; Seitz & Schneider 1997). The initial mass re- 
construction can be carried out by setting k — 0. Now, since the improved information on 
k is available, we can estimate the amount of correction for the original reduced shears g, 
and the updated 7 can be resubmitted into the mass reconstruction. This k — > 7 — > k cycle 
converges typically after a few iterations. The other more fundamental problem is that the 
shear in equation 23 is invariant under k — > Xk + (1 — A) transformation. This so-called 
sheet mass degeneracy cannot be lifted unless external constraint is provided. Furthermore, 
without the proper knowledge of this rescaling, we cannot achieve the previous nonlinear 
mass reconstruction because the validity of k for the use of updating shear 7 is not guaran- 
teed. In the next section we demonstrate that this rescaling of the surface density k map 
is in reality possible with the help of the parameterized cluster profile. The parameterised 
mass modeling is safe from this sheet-mass degeneracy though the accuracy is sometimes 
compromised due to the inadequacy of the assumption of a particular mass profile. We will 
also show that the direct mass estimates from this "rescaled" k map are in good agreement 
with the results from the conventionally favored aperture densitometry (Fahlman et al. 1994; 
Clowe et al. 1998). 



7. MASS ESTIMATES 

In the current section we present the mass estimates of CL 0152-1357 through vari- 
ous routes. First, we discuss the results from parameterized profile fitting. Especially for 
high-redshift clusters, which have a significant substructure, one does not expect that pa- 
rameterized models can optimally describe complex profiles. Nevertheless, parameterized 
model fitting is an invaluable procedure not only because it can easily provide a reasonable 
first-guess of the mass profile of the cluster, but also because the results can be used to 
estimate the feasible mean surface mass density of the annulus far from the center, thus en- 
abling one to constrain R of this region in subsequent mass estimation. Second, we consider 
aperture mass densitometry which has been a preferred choice in many cluster weak lensing 
studies because its implementation is straightforward and safe from artifacts arising in most 
reconstruction algorithms. Finally, we attempt to measure the mass of the cluster directly 
from the mass reconstruction map presented in the previous section. We show that, after 
the consideration of proper rescaling k — > Xk + (1 — A) and non-linearity g — 7/(1 — k), 
the mass estimate obtained from the mass map is very close to the results from aperture 
densitometry. 
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7.1. 



Parameterised Mass Profile 



7.1.1. Singular Isothermal Sphere and Ellipsoid 



For a Singular Isothermal Sphere (SIS), the dimensionless surface density k is given as 



where Be is the Einstein radius. It is also easy to show that the simple relation |7| = k exists 
for SIS. After iteratively transforming the observed (reduced) tangential shears into the true 
tangential shears 7 = (1 — n)g, we fit a SIS profile to those tangential shears in the annulus 
from 50" to 160". From the typical \ 2 minimization we find 9 E = 6.64 ± 0.82". In addition, 
we test if the ellipticity of the cluster can be detected by fitting a Singular Isothermal 
Ellipsoid (SIE) (Kormann, Schneider, & Bartelmann 1994). While a single parameter can 
characterize SIS, SIE requires 3 parameters which describe Einstein radius 9e, axial ratio 
/, and orientation angle a. Following the notation of King & Schneider (2001), the surface 
mass density and shear are related as follows: 



where / and a are axis ratio and orientation angle of a cluster, respectively. The correspond- 
ing shears are 



For SIE profile fitting we use the smoothed distortion field in the same annulus as in SIS 
fitting. The best-fit parameters are 6 E = 7.19 ± 0.72", / = 0.36 ± 0.10, and a = 17 A ± 6.7° 
(from the vertical axis). The orientation angle is consistent with the distribution of the 
cluster galaxies as well as the mass reconstruction. The axial ratio / indicates that the 
overall mass distribution is highly elongated and the azimuthal variation of the surface mass 
density is still detectable even at large radii. We determine the velocity dispersions from the 
estimated Einstein radii of SIS and SIE to be of JS = 903+^ km s" 1 and a% IE = 940+^ km 
s _1 , respectively. These results are in good agreement with the direct measurements from the 
redshift survey of the cluster (R. Demarco et al. 2004, in preparation). They measured the 
velocity dispersion from their spectroscopic data within ~ 50" aperture radii and obtained 
o~ v = 919 ± 168 km s _1 and o~ v = 737 ± 126 km s _1 in rest-frame for the northern and 
southern clumps, respectively. It is also possible to compare the velocity dispersions with 
the cluster temperature estimates from the X-ray and Sunyaev-Zeldovich analyses. Despite 
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71 = — kcos(20) 

72 = — /tsin(20). 



(27) 
(28) 
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the fact that the a — T relation of the galaxy clusters is in general very scattered along 
the theoretical prediction a v oc T 1 / 2 line, the estimates of 5.9 +4 ' 4 keV (Ebeling et al. 2000), 
6.5+{;J; keV (Delia Ceca, Maccacaro, Rosati, k Braito 2000), 8.5±f;g keV (Joy et al. 2001), 
or 5.5+Q g keV (Maughan et al. 2003) are consistent with these velocity dispersions. For 
example, if we adopt the empirical relation (cr^/krn s~ x ) = io 2 - 57±ai3 (/cT/keV)°' 59±0 ' 14 (Wu 
et al. 1998), we get T ~ 4. 51^ keV for the SIS velocity dispersion. 

7.1.2. NFW Profile Fitting 
The NFW density profile (Navarro, Frenk, & White 1997) is defined as: 



o( r ) - 5 J*M (29) 

P[) (r/r s )(l + r/r s r 1 J 



where 8 C is the halo overdensity expressed in terms of the concentration parameter c as 

200 c 3 



3 ln(l + c) - c/(l + c) ' 



(30) 



p c (z) is the critical density of the Universe at the redshift of the cluster, and r s is the scale 
radius of the profile. The relation between mean surface density and gravitational shear 
is much more complicated in the NFW profile. The useful mathematical formalisms for 
gravitational lensing are worked out by Bartelmann (1996), Wright & Brainerd (2000), and 
King & Schneider (2001). 

Though many parameters seem to be involved in the characterization of the NFW 
profile, only two parameters are independent. From the similar x 2 minimization as in SIS 
and SIE fitting, we estimate the concentration parameter c and the scale radius r s to be 
3.7 ± 0.5 and 40 ± 6" (309 ± 45 kpc), respectively. A virial radius is defined as a radius 
where the mean interior density drops to 200 times the critical density of the Universe at 
the redshift of the cluster. From the simple relation r 2 oo — cr s , the virial radius is estimated 
to be r 200 = 150 ± 30" (1.14 ± 0.23 Mpc). 



7.2. Aperture Mass Densitometry 



When one's interest is to find a total mass within some given aperture radius r, the 
following ((r) statistics provide a useful measure of lower limits on the mean surface mass 
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density inside r: 



( c (n,r 2 ,r. 




dr + 




<7t) 



dr, 



(31) 



max 



1 — r; 



r 



where {•jt) is an average of tangential shears defined in equation 19, 77 is the aperture radius, 
and r 2 and r max are the inner- and the outer radii of the annulus. In the weak lensing 
regime where the relation between the measured ellipticity and the true shear becomes 
linear, Cc(?"i, r 2 , r max ) can be directly computed from the tangential shear to estimate R(r < 
77) — R(r 2 < r < r max ) where re(r 2 < r < r max ) is an average mass density in the annulus. 
In principle, if one choose r 2 and r max far enough from the cluster center to make the 
contribution vanishingly small, the ((r) above approaches the genuine average surface mass 
density within the radius 77. In the present study, we used r 2 = 140" and r max = 160" 
in order to keep the entire annulus within the observed field. The mean surface density in 
this annulus is expected to be low, but it still contributes to £(r). From the result of the 
SIS fit, we estimate the dimensionless mean surface density of this region R(r 2 < r < r max ) 
to be 0.023 ± 0.003. Figure 16 shows the mean surface density inside given radii after the 
contribution from the annulus R(r 2 < r < r max ) is added. At small radii, re is overestimated 
because of the reasons discussed in §6.2 and §6.4. The dashed line represents the mean surface 
density when this correction is applied. The difference amounts to ~ 11% at r ~ 20". 



In 6.4 we discussed the general difficulties in translating the reconstructed convergence 
map into the mass density in physical units. In order to lift the re — > Are + 1 — A degeneracy, 
we must be able to constrain re at least for a limited region of the field. This is not impossible, 
however, because in §7.2 we were able to estimate the mean surface density R(r 2 < r < r max ) 
in the control annulus from the parameterized mass models. Therefore, we can compare 
this with the value measured in the same annulus of the mass map. The transformation 
parameter A becomes no longer arbitrary and can be determined from the relation 



where Rl is a mean surface density of the annulus from the parameterised models and R is the 
same quantity measured from the mass map. Then, we can apply Are + 1 — A transformation 
to the entire region of the mass map. Since this mass map is properly scaled, we can use 
it to update the shear and feed this corrected shear back to the mass reconstruction. The 
procedure is iterated until convergence is reached. In this study, no more than 4 iterations 
were needed. 



7.3. Rescaling of the Mass Reconstruction Map 



A 



re- 1 



(32) 
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We compare the cumulative projected mass profile from this rescaled mass map with the 
result from the aperture densitometry in Figure 17. The excellent agreement between these 
two profiles demonstrates that at least in an azimuthally averaged sense the mass estimate 
from the mass map is consistent with the result from the aperture densitometry. In §8, we will 
also determine the mass of the sub-clumps through both of these routes and show that the 
consistency can be generalized. As far as the rescaling is appropriately calibrated, the use of 
the rescaled mass map in probing the mass distribution of the cluster has obvious advantages 
over aperture densitometry. The aperture densitometry always requires the control annulus 
to be set up around target apertures, which is sometimes hindered if the annulus cannot form 
a complete circle. Besides, to estimate the mean surface mass density inside the annulus, one 
has to fit a particular parameterised mass model. If the discrepancy between the assumed 
and the actual cluster mass profile is not small (e.g. due to the substructure), the procedure 
always introduces additional uncertainties. On the contrary, the direct use of the rescaled 
mass map does not suffer from these obstacles, and this method can become particularly 
useful when mass inside some arbitrary boundary needs to be estimated. 



7.4. Comparison Between Parameterised and Parameter-Free Methods 

The mass profiles from the SIS, NFW, and aperture densitometry are compared in 
Figure 18. We omit the SIE fitting result because it overlaps the SIS profile very closely. 
Obviously, the actual cluster mass profile is best approximated by the NFW profile (solid). 
If the SIS (dotted) is assumed instead, the total projected mass is overestimated by ~ 20% 
at r ~ 131" (~1 Mpc). Considering the apparent filamentary substructure of CL 0152- 
1357 delineated by either the light or the mass distribution, the excellent representation 
of the azimuthal mass distribution of the cluster by the NFW profile is rather remarkable. 
We summarize the mass estimates inside 1 Mpc radius aperture in Table 1 for various 
combinations of cosmological parameters and methods. 



7.5. Mass-to-light Ratio 

We present the mass-to-light ratio profile of CL 0152-1357 in Figure 19. The cumulative 
mass-to-light ratio M(< r)/ L B &{< r) (open circle and dashed) of the cluster rapidly rises to 
its maximum at r ~ 35" and then decreases rather monotonically. The decrease of the profile 
looks more pronounced in the differential mass-to-light ratio 5M{r) / 5Lbq{t) (dotted). It is 
verified that the profile when the blue cluster galaxies are excluded does not significantly 
change though the slope is slightly reduced. The small mass-to-light ratio near the cluster 
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center seems to originate from the luminosity segregation of the brightest cluster galaxies. 
Carlberg, Yee, & Ellingson (1997) studied the average mass-to-light profiles of 14 galaxy 
clusters from the virial mass estimator. The resulting mass-to-light ratio averaged over 
substructure and asymmetries is high in the inner regions and gradually decreases until 
it starts to flatten at r ~ 0.7r2oo- The average mass-to-light profile obtained from the 
kinematics and distribution of 3056 galaxies in 59 nearby clusters in the ESO Nearby Abell 
Cluster Survey also shows a similar trend of a rapid rise followed by a gradual decrease up 
tor ~ 0.7r 20 o- Does the upturn of the M/L profile of CL 0152-1357 at r ~ 115" correspond 
to the beginning of the plateau observed in those works? Assuming the feature is real and 
the empirical relation r 2 oo — r/0.7 holds, the virial radius of the cluster can be evaluated to 
be r 2 oo ~ 164" (~ 1.2 Mpc). This value is surprisingly close to the independent estimation 
from NFW fitting in §7.1.2. 

The average M/L ratio of the cluster within 1 Mpc aperture radius is estimated to be 
95 ± 8 Mq/Lbq- It is of interest to compare this value with the result for MS 1054-03 at a 
very similar redshift of z ~ 0.83. In Qm = 0.3 and Q\ = 0.7 cosmology HFK00 quoted the 
M/L ratio of 127±10h 71 M Q /L BQ , which is higher than that of CL 0152-1357 by ~ 34%. The 
result is consistent with the still arguable, but popular belief that the cluster M/L increases 
with richness. The M/L profile of MS 1054-03 by HFK00 shows the gradual increase of the 
M/L ratio out to the field limit. If the M/L profile of MS 1054-03 is assumed to conform to 
the aforementioned average M/L profile at large radii, we may suggest that the aperture of 
1 Mpc in MS 1054-03 field encompass only the inner region where the M/L is still high. 

The logarithmic luminosity evolution, In [Mj Lb) oc (—1.06 ± 0.09)z, is derived by van 
Dokkum & Stanford (2003) from massive cluster galaxies at 0.02 < z < 1.27. At z — 0.84, 
the relation predicts ~ 41% reduction in B-band luminosity and the M/L ratio of CL 0152- 
1357 is modified to be ~ 232M Q /L B0 , which is similar to the results for other clusters (e.g., 
Carlberg, Yee, & Ellingson 1997). 

8. Substructure of CL 0152-1357 

8.1. Mass Estimates of Individual Mass Clumps 

Due to the high number density of background galaxies whose shapes are reliably mea- 
surable, the reconstructed mass map reveals the cluster substructure in great detail. We 
overplot the rescaled mass reconstruction on the negative gray image of CL 0152-1357 in 
Figure 20. We identified 9 mass clumps whose significance is above ~ 3<x and galaxy coun- 
terparts are apparent. The significance for each mass pixel is computed by the use of the 
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RMS mass map (Figure 21), which is constructed from bootstrap 5000 realizations of mass 
reconstruction. These 5000 mass maps are also used in §8.2 to examine the uncertainties 
of the mass peak centroids. The mass of these clumps within 20" aperture (~ 150 kpc) 
are computed via the direct use of the reconstruction map as well as the examination of 
£(r) statistics (Table 2). They are in good agreement with each other with overlapping 
uncertainties. 



8.2. Comparison with Other Studies 

Though the Einstein IPC first detected X-ray emission from CL 0152-1357, its signifi- 
cance was not properly recognized because of the complex morphology of the emission. The 
cluster was rediscovered in the WARPS (Scharf et al. 1997), the RDCS (Rosati, della Ceca, 
Norman, & Giacconi 1998) and the SHARC (Romer et al. 2000) survey. Ebeling et al. (2000) 
analyzed the X-ray observation of CL 0152-1357 from the WARPS survey and showed that 
the X-ray morphology of the cluster is suggestive of very complex substructure that can be 
also traced by cluster galaxy concentrations. The higher resolution of Chandra extended 
the work by Ebeling et al. (2000) and revealed two prominent X-ray peaks (Maughan et 
al. 2003). In order to verify their results and also enable a direct comparison of the X-ray 
morphology with our weak lensing mass distribution, we reanalyzed the archival Chandra 
observations. After adaptively smoothing the X-ray image via "csmooth" , which is part of 
the Chandra Interactive Analysis of Observations Software (CIAO), we obtained the X-ray 
contour map of the cluster which we overlaid on the ACS image (Figure 22). The astrometric 
accuracy of the ACS image with respect to the Guide Star Catalog 2 (GSC-2) is ~ 0.05". 
Considering ~ 0.1" errors on the GSC-2 itself and also the absolute astrometric accuracy of 
the Chandra observation being ~ 1", we expect the alignment between the two images are 
fairly precise, and the excellent agreements of the locations of X-ray point sources with those 
of the optical counter parts confirm this fact. The two main diffuse X-ray peaks, though 
slightly off-center from two concentrations of spectroscopically confirmed cluster members, 
seem to correspond to the mass clumps C and F in our weak lensing map (Figure 20). The 
mass clump A, which is associated with the z = 0.846 ±0.003 group of galaxies (R. Demarco 
et al. 2004, in preparation), is also detected as extended low surface brightness region in 
the Chandra observation though this location is somewhat remote from that of the galaxy 
concentration. In addition, one of the point-like X-ray sources appears to be associated 
with the mass clump E. Maughan et al. (2003) examined the X-ray excess emission from 
the residual image which is constructed by subtracting the best-fit model. They discovered 
that the regions of excess emission lie midway between the two major clumps, stretched 
perpendicular to the merging direction. They also found that the X-ray temperature of this 
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region is relatively higher. Though they were not conclusive because of the low significance 
of the feature, the presence of a shock front was suggested. It is remarkable that the feature 
is also seen in the weak lensing mass map as clumps D and E. However, the cause of the 
mass concentration in this region cannot be exclusively ascribed to a shock which originates 
from the merger of two major subclusters. The cluster light distribution suggests that there 
may also be cluster galaxies associated with these clumps. 

It is interesting to examine the displacements of peaks between the weak lensing mass 
map and the X-ray flux contours. We present overplots for three different combinations: X- 
ray/optical (Figure 23), mass/optical (Figure 24), and mass/X-ray (Figure 25). In general, 
weak lensing mass distribution better traces the light distribution of the cluster as far as 
the overall morphology and coincidence of clump locations are concerned. As Maughan et 
al. (2003) observed, the location of the southern X-ray peak is offset by ~ 5" from the 
galaxies and displaced away from the direction of the merger. They proposed that the 
displacement is due to the fact that relatively collisionless galaxies are moving ahead of the 
viscous ICM whose distribution is delineated by the X-rays. We observe that the similar 
trend for the northern X-ray peak is also noticeable but with smaller displacement. The 
offset is more obvious when the location of the northern X-ray peak is compared with the 
smoothed luminosity center (not with the brightest galaxies; Figure 23). By examining 
the significance of the event counts around the northern X-ray peak, we verified that the 
offset is not likely to be caused by artifacts of adaptive smoothing. Similar, but more 
pronounced centroid offsets were detected in the Chandra study of the merging cluster IE 
0657-56 (Markevitch et al. 2002). The X-ray image shows that the X-ray centroids of two 
clumps are conspicuously displaced from the cluster galaxy concentrations, suggesting the 
ram-pressure stripping of the gas components of the cluster. If the dark matter is indeed 
collionless and is the dominant contributor to the total cluster mass, the location of the mass 
peaks must be separated from the X-ray centroids. Clowe, Gonzalez, & Markevitch (2004) 
performed a weak lensing analysis of the same cluster and showed that the reconstructed 
mass peaks are displaced from the X-ray halos while in good spatial agreements with the 
galaxy concentrations. 

In our mass reconstruction of CL 0152-1357, we find that the two mass clumps detected 
in the weak lensing map are also shifted toward the suspected merging direction with respect 
to the luminosity as well as the X-ray distribution (Figure 24 and 25). Centroids of weak 
lensing mass clumps are in general affected by shape noise, shear strength, reconstruction 
algorithm artifacts, spatial number density variations of source galaxies, etc. Though most 
of these factors are present to a varying degree in our mass reconstruction, we expect the 
uncertainties of these two centroids are relatively low compared to those of outside the main 
body because of their high significance in the reconstructed mass map. Furthermore, it 
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is understood that the high number density of source galaxies (~ 175 objects arcmin~ 2 ) 
increases the overall stability of the centroids of our weak lensing mass peaks. One of the 
useful tests to quantify the centroid uncertainty is to examine how much the locations of the 
mass peaks change when the shears are perturbed slightly. Motivated by the experiment of 
Clowe, Gonzalez, & Markevitch (2004), we measured the centroids of the two mass peaks in 
5000 runs of mass reconstruction obtained in § 8.1 by bootstrap resampling. The resulting 
distribution of the two mass peaks are presented in Figure 26. It shows that both of the 
two light peaks are outside the 99% circle. The RMS distance of the northen (clump C) 
and southern (clump F) peaks are ~ 0.4" and ~ 1.0", respectively. The smaller RMS 
scatter of the northern peak is consistent with the higher significance in the reconstructed 
mass map. The presence of these preferential shifts in the detected clumps may extend the 
argument above and further support the collisionless nature of CDM. Within the paradigm 
of the hierarchical structure formation, the shifts may indicate that the CDM which initially 
created the deep potential well for the formation of the subclusters is moving even faster than 
galaxies which are less subject to ram pressure than the ICM, but not entirely collionless as 
the CDM. 

We observe that there lie four bright foreground (z ~ 0.63) galaxies to the south of 
the BCGs. If they are massive enough to perturb the distortion of background galaxies, 
the centroid of the northern clump can be affected and appear to be further shifted toward 
the merging direction. However, though we cannot completely rule out this possibility, we 
suspect that their contribution to the centroid shift is not so substantive as to cause the 
distinct offset between the X-ray and dark matter contours (Figure 25). More study of these 
questions regarding the centroid offsets among galaxies, X-ray emission, and weak lensing 
mass will be conducted when the weak lensing analysis of another supposedly merging cluster, 
MS 1054-03 at z ~ 0.83, with ACS observations (M. Jee et al. 2004, in preparation) becomes 
available. 

We compare the mass estimates of the cluster with those from the Chandra X-ray 
spectral analysis (Maughan et al. 2003) and the Sunyaev-Zeldovich Effect (SZE) work (Joy 
et al. 2001) by treating the cluster as a whole to simplify the comparison. Joy et al. (2001) 
infer that the total mass within a radius of 65" in fl M = 0.3, fl\ = 0.7, and H Q = 100 
kms _1 Mpc _1 universe is (2.1 ± 0.7) x 10 14 M Q from their SZE measurement of the cluster 
temperature. The weak lensing mass under the same cosmological parameters yields (1.9 ± 
0.2) x 10 14 M Q . From the Chandra X-ray analysis Maughan et al. (2003) quotes an estimate 
of 2Atoi x 1O 14 M within ~ 50" radius aperture in Vt M = 0.3, fi A = 0.7, and H = 70 
kms _1 Mpc~ 1 universe. Our conversion of the weak lensing mass under the same geometry 
is estimated to be (2.1 ± 0.3) x 1O 14 M . We note that the weak lensing mass estimate is 
lower by ~ 9% in both comparisons though the statistical significance of the difference is 
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low. In the temperature-based measurement of the cluster mass, the dynamical equilibrium 
of the intracluster gas with the underlying dark matter as well as the spherical symmetry of 
the matter distribution must be assumed. However, the on-going merger obviously violates 
these hypotheses and especially the temperature rise caused by the shock may result in the 
overestimation of the cluster mass. 

The X-ray mass estimate of (1.1 ±0.2) x 10 15 M Q within 1.4Mpc radius (extrapolated to 
the virial radius under the assumption of the isothermal sphere) by Maughan et al. (2003) 
is roughly a factor of two higher than the result from the current paper even considering the 
aforesaid risk of overestimation as well as the dissimilar geometry. We attribute this rather 
large difference to the two following reasons. First, as discussed in §7.4, the mass profile 
of CL 0152-1357 is better described by the NFW profile. The SIS modeling of the cluster 
profile gives substantially higher total mass than the NFW representation does at large radii 
(leading to ~ 37% increase at ~ 1.4 Mpc). Second, the mass of the southern X-ray peak 
from this weak lensing analysis is much lower than the northern one (less than 50 percent 
of the northern X-ray peak within 20" aperture radius, see Table 3) whereas Maughan et al. 
(2003) estimates that the two peaks are of similar mass. Therefore, the total virial mass of 
CL 0152-1357 computed by the superposition of two comparable SIS clumps is likely to be 
much higher than the result from the current analysis. 

Huo et al. (2004) presented the first weak lensing analysis of CL 0152-1357 using Keck R 
band observations. The projected cluster mass of ~ 1O 15 M within r = lh^ Mpc can be read 
off Figure 10 of their paper. The transformation of our mass estimate in their cosmological 
parameters gives (5.09±0.46) x 1O 14 M . Despite the somewhat large discrepancy in mass, we 
do not further analyze the difference because the intermediate procedures of the weak lensing 
analysis (e.g., PSF corrections, tangential shears, reconstruction maps, redshift distribution 
of background galaxies, etc.) are not illustrated in their work. 

9. SUMMARY AND CONCLUSIONS 

We have presented our weak lensing analysis of the X-ray luminous cluster CL 0152- 
1357 at z ~ 0.84 using ACS observations. The superb resolution and sensitivity of ACS 
provides high quality images of weakly distorted, faint background galaxies in unprecedented 
depth. The resulting high number density of source galaxies enables us to restore the cluster 
mass distribution in unparalleled detail when the instrument artifact is properly accounted 
for. The complicated shape and variation of the PSF is precisely modeled by exhaustive 
investigation of the 47 Tuc stars, and the derived PSF is matched to the isolated good signal- 
to-ratio stars in CL 0152-1357 field after a slight fine-tuning of the ellipticity. Rounding kernel 
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test shows that the final PSF obtained in this way nicely describes the observed PSF pattern 
in CL 0152-1357 field. We use the publicly available GOODS and UDF images to infer the 
fraction of cluster galaxy contamination in source galaxies, and the redshift distribution of 
background galaxies is estimated by the use of the photometric redshift catalog of the UDF. 

We determine the cluster mass via three different approaches: parameterized profile 
fitting, mass reconstruction, and aperture mass densitometry. Among these approaches, 
the second method of mass estimation from the reconstruction map is unconventional. The 
direct use of the weak lensing mass map has been discouraged primarily because there exists 
a sheet-mass degeneracy. However, we show that, after proper rescaling is considered, the 
method yields very consistent results with the measurements from the often favored aperture 
mass densitometry. 

Our weak lensing mass estimates at small radii (r < 65") are consistent with the results 
from the X-ray emission and the Sunyaev-Zeldovich effect. We show that the deviation of 
the mass profile from the SIS profile increases at larger radii. The overall mass profile of the 
cluster can be well described by the NFW profile with a scale radius of r s = 309 ±45 kpc and a 
concentration parameter of c = 3.7±0.5. We estimate the total projected cluster mass within 
1 Mpc aperture to be (4.92 ± O.44)xlO 14 M from the aperture mass densitometry. The total 
luminosity of the cluster is calculated by combining the spectroscopic and the red sequence 
catalogs. When 2775 is transformed to the rest frame B, the total luminosity (r < 1 Mpc) after 
accounting for the blue and the faint population is determined to be Lbq = 5.2 x 10 12 Lb & . 
We find that the M/L ratio within 1 Mpc is 95± 8 M q /Lbq- Considering the luminosity 
evolution at z ~ 0.84, this M/L ratio corresponds to ~ 232M Q / Lbq at local universe. 

Our ACS weak lensing reveals very interesting substructure of the cluster in detail. The 
vertically elongated cluster main body is clearly seen in both light and mass distributions with 
strong spatial correlations between light and mass clumps. Besides, we identify 4 scattered 
mass clumps outside the main body with locations of cluster galaxy concentrations. More 
stimulating interpretation is made when the mass reconstruction is compared with the X- 
ray morphology from Chandra observations. In order to examine the spatial correlations 
between the two analyses, we reprocess the Chandra archival data and overlay the X-ray 
contours with those of optical light and mass maps. We observe that the two diffuse X-ray 
clumps are in spatial agreement with cluster galaxy concentrations, but are displaced away 
from the assumed merging direction. The displacement of the southern peak was originally 
noticed by Maughan et al. (2003) and they suggested that the ICM is lagging behind the 
cluster galaxies due to the ram pressure. The comparision of the X-ray emission with our 
mass reconstruction strengthens this merger hypothesis because both the cluster galaxies and 
mass clumps seem to lead the X-ray peaks. Furthermore, we remark on the displacements of 
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the mass clumps relative to the light concentrations. It appears in the mass-light comparison 
map the two major mass clumps are slightly shifted (~ 3") toward the merging direction. 
The existence of these preferential shifts might suggest that the collisionless dark matters 
are moving ahead of cluster galaxies. 

Though the weak lensing survey data from today's extensive dedication of many large 
aperture ground-based telescopes, primarily targeted for the cosmic shear detection, surpass 
those of ACS in data volume, our weak lensing analysis of CL 0152-1357 demonstrates that 
ACS is exclusively advantageous for weak lensing studies of high-redshift clusters which 
require only moderately large field of view, but extremely high resolution of the instrument. 
The ACS GTO cluster survey program encompasses many high-reshift clusters of great 
interest and the weak lensing investigation of these clusters will provide many illuminating 
clues to the formation and evolution of galaxy clusters in the near future. 

ACS was developed under NASA contract NAS 5-32865, and this research was sup- 
ported by NASA grant NAG5-7697. We are grateful for an equipment grant from the Sun 
Microsystems, Inc. 
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Fig. 1. — Color composite of CL 0152-1357 center. The image shows the central ~ 40" x 40" 
section of the entire ~ 350" x 350" field. The image is created using the FITSCUT (McCann 
2004). 
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Fig. 2. — Number counts of galaxies in the r 62 5 (solid), i 775 (dotted), and z 850 (dashed) 
passbands. The detection is complete down to ~ 27.5 mag in all filters. 
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Fig. 3. — Statistics of galaxies whose shapes are reliably measured. Solid line (top panel) rep- 
resents the number of objects detected by SExtractor and dashed line indicates the number 
of objects for which the iteration converges. The fraction of measureable galaxies decreases 
substantially after ~ 29 mag (see bottom panel). 



-35 - 



4000 



3000 

tn 

JL 2000 

>- 

1000 







■I'" I 1 " I 1 



fnV 



, , " ' "I H 



- - Sua i „ - ' «t.\ . '1 



1000 2000 -5000 4000 
X (plxeis) 



4000 



.5000 

m 

J, 2000 



1000 ■ 







1 ' 1 ■ 1 


i ' ■ 




! 


■ 

< 








■ 1 




!. " _ 1 ' 

I ; 

, 1 1 ' 1 ' 1 - 
r i i ' ■ ™ 



1000 2000 .5000 4000 
X (pixeis) 



Fig. 4. — PSF variation of WFC sampled from external field observation of 47 Tuc. The 
length and direction of the whisker are proportional to the magnitude and orientation angle 
of ellipticity, respectively, (a) PSF pattern observed on 2002 October 3. Most stars are 
elongated from lower-left to upper-right.; (b) The same field is taken on 2002 October 24. 
The PSF elongation is nearly perpendicular to the pattern in (a). 
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Fig. 5. — PSF pattern in r 62 5 passband of CL 0152-1357 field, (a) Uncorrected ellipticities 
measured from stars, (5) = 0.114 ± 0.040; (b) Ellipticities after removal of PSF anisotropy 
(5) = 0.023 ± 0.017. Correction is made by reconvolving the image with position-dependent 
rounding kernels. The circled whiskers at center illustrate 10% ellipticity. 
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Fig. 6. — PSF anisotropy removal on CL 0152-1357 field. Diamonds represent the initial 
ellipticities of stars and plus symbols the corrected ellipticities after rounding kernel convo- 
lution. The improvement is remarkable in both amplitude (size of scatter) and anisotropy 
(centroid). 
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Fig. 7.— Color-magnitude (CM) plot of CL 0152-1357. The tight CM relation of the early- 
type cluster galaxies is present at (r 62 5 — z$5o) ~ 1.9. Bright stars are not removed. 
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Fig. 8. — Luminosity map of CL 0152-1357. For galaxies brighter than z 850 = 22, spectro- 
scopically confirmed members are selected whereas down to Zs^o — 25 member selection is 
based on the color-magnitude relation of early-type galaxies. The luminosity center does not 
exactly lie on the BCGs, but is slightly shifted to the south by ~ 2" because of the presence 
of other bright cluster galaxies. 
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Fig. 9. — Cumulative luminosity profile. The i 775 band magnitudes are transformed to 
the rest-frame B band absolute magnitudes after corrections are made for the galactic dust 
extinction and the distance modulus. 
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Fig. 10. — Shear recovery test. The simulation is performed by artificially shearing the r^s 
image of CL 0152-1357 by 7 = 0.05 and measuring the difference in ellipticities of individual 
galaxies. The large uncertainty for bright galaxies is due to the reduced number of objects. 
It is verified that down to ~ 29.5 mag the input shear is successfully recovered with less than 
£7 ~ 0.01 deviation. 
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Fig. 11. — Tangential shear versus radius from the cluster center. (gj<) for three different 
samples. The tangential shear from bright sample (z 850 < 26; top) is weak in comparison 
with those from faint (z 850 < 28.5; middle) or faintest (z 850 < 30; bottom) sample. We 
also display best-fit singular isothermal sphere models (dashed) with corresponding Einstein 
Radius (0 E )- The inner region (r < 50") is excluded from the fit. 
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Fig. 12. — Cross shear test for gravitational lensing. Tangential shear is evaluated after 
source galaxies are rotated by 45° (top). The absence of the signal verifies that the tangential 
shear in Figure 11 is due to the gravitational lensing. The comparison of this result with the 
shuffle test (bottom) where galaxy ellipticities are randomly shuffled while the positions are 
fixed demonstrates that the cross shear scatters are consistent with the randomization test 
results. 
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Fig. 13. — (a) Contamination by cluster galaxies. The magnitude distribution of the source 
sample (solid) is compared with those of the GOODS (dotted) and the UDF (dashed). 
Because vqoq filter is used in both surveys instead of r 62 5, we transformed (t> 60 6 — -2850 ) color 
to ( r 625 — z 85o) color (M. Sirianni et al. 2004, in preparation) to apply a consistent selection 
criterion. The fraction of cluster galaxies in the source sample can be estimated from the 
magnitude distribution of the GOODS galaxies down to z 850 ~ 26. At z 850 > 26, the 
incompleteness of the GOODS catalog is clear. In order to estimate the contamination in 
this magnitude range, we degrade the UDF images to mimic the S/N ratio of our cluster 
observation and detect source objects on these simulated images (open circle), (b) Estimation 
of (/3) as a function of magnitude. The contamination by cluster members are not included 
yet in this plot. 
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Fig. 14. — Shear field of CL 0152-1357. Galaxy ellipticities measured from different passband 
images are combined and smoothed with a Gaussian kernel (FWHM ~ 20"). 
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Fig. 15. — The convergence k map reconstructed from the shear grid. The shear invariant 
transformation k — > A/t + (1 — A) is not applied yet. The positions of distinct mass clumps 
are highly correlated with the luminosity peaks found in Figure 8. 
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Fig. 16. — Projected mean surface density within apertures in units of critical density S cr . 
«(140" < r < 160") = 0.023 ±0.004 is added to the ((r) statistic. The dotted line represents 
the same mean surface density when the (7 = 7/(1 — k) correction is included. This factor is 
important at the cluster core where the k C 1 assumption breaks down. 
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Fig. 17. — Enclosed projected mass profile. The mass estimates derived from the mass 
reconstruction are compared with those from the aperture densitometry (open circles with 
error bars). The first, second, and final (4th) iterations of rescaling of the mass map are 
represented by solid, dotted, and dashed lines, respectively. 
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Fig. 18. — Comparison of the parameter-free mass profile with the results from SIS and NFW 
fitting. The parameter-free cluster mass profile (open circle with error bars) is excellently 
described by the NFW profile with a scale radius of r s = 40 ± 6" (309 ± 45 kpc) and a 
concentration parameter of c = 3.7 ± 0.5. The SIS profile (dashed) overestimates the cluster 
mass by ~ 20% at r ~ 131" (~1 Mpc). 




Fig. 19. — Mass-to-light ratio profile of CL 0152-1357. Open circles and the dashed line 
represent the M/L ratio based on the aperture mass densitometry and the reconstructed 
mass map, respectively. The dotted line shows the differential M/L ratio. The uncertainties 
shown do not include the photometric errors for individual galaxies. The cumulative M/L 
ratio peaks at ~ 35" and then decreases rather monotonically, reaching ~ 95M / Lbq at 1 
Mpc (~ 131") from the cluster center. This trend is more distinct in the differential M/L 
profile. The light profile (Figure 9) is smoothed to reduce the severe scatters caused by the 
discrete galaxy positions. 
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Fig. 20. — Overlay of the mass map on the negative gray image. We only displayed contours 
for k > 0.08 which corresponds to ~ 3o\ The mass estimates of annotated clumps are 
summarized in Table 2. 
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Fig. 21. — RMS map of the rescaled mass map in Figure 20. The rms is computed by 5000 
runs of mass reconstruction. For each run, we randomly resampled background galaxies 
from our source catalog. Note that at the field boundaries (especially at the four corners) 
the uncertainties of mass pixels are higher than in the inner region. 
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Fig. 22. — Overlay of the smoothed X-ray map on the ACS detection image. The relative 
astrometric accuracy between the Chandra and the ACS images is ~ 1.4" (see text). The X- 
ray flux map is generated after adaptively smoothing the raw X-ray image using "csmooth" 
which is part of the Chandra Interactive Analysis of Observations Software (CIAO). 
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Fig. 23. — Overlay of the smoothed X-ray contour (white solid) on the 2775 luminosity 
distribution (color-coded). The X-ray peaks of two subclusters seem to be displaced away 
from the assumed merging direction. 
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Fig. 24. — Overlay of the mass contour (white solid) on the 2775 luminosity distribution 
(color-coded). The locations of mass clumps are in contrast to the X-ray peaks in Figure 23, 
shifted toward the cluster center. 
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Fig. 25. — Overlay of the weak lensing mass map (white solid) on the smoothed x-ray 
background (color-coded). The offsets between mass clumps and X-ray peaks are distinct. 



-4-2 7 4 

X (crcsec) 



-6 -A -2024 I 
X (crcsec) 



Fig. 26. — Distribution of centroids in 5000 mass reconstructions. The luminosity peaks are 
marked by plus (+) symbol. In both figures, the X-ray peaks are outside the plot range. 
The circle represents the area where 99% of the data is enclosed. No rejection is made, (a) 
Centroid distribution of the northern peak (clump C). The RMS is ~ 0.4"; (b) Centroid 
distribution of the southern peak (clump F). The uncertainty increases to ~ 1.0". 
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Table I. MASS PROPERTIES OF CL 0152-1357 (r < 1 Mpc) 





n A 


H 


M S IS 




M <(r) 




M ((r) /L 






(km/s/Mpc) 


(1O 14 M ) 


(1O 14 M ) 


(1O 14 M ) 


(1O 14 M ) 


(M q /L B q) 


0.27 


0.73 


71 


6.07 ± 0.75 


4.92 ± 0.71 


4.92 ± 0.44 


4.74 ± 0.12 


95 ± 8 


1.0 





50 


7.31 ± 0.90 


6.05 ± 0.88 


5.94 ± 0.54 


5.86 ± 0.12 


108±10 



Table 2. PROPERTIES OF SUBCLUMPS IN CL 0152-1357 



Subclump 


C(< 20") 


K(r2 <r< Tmax) 


r-2 


'r'max 


M C(r) 


Mm.ap 


Mmap / L 








(") 


(") 


(1O 13 M ) 


(1O 13 M ) 


(M @ /L B q) 


A 


0.064 ± 0.026 


0.019 ± 0.013 


40 


60 


2.2 ± 0.8 


2.1 ± 0.3 


73 ± 12 


B 


0.120 ± 0.027 


0.022 ± 0.006 


80 


100 


3.7 ± 0.7 


3.4 ± 0.2 


80 ± 6 


C 


0.232 ± 0.032 


0.023 ± 0.003 


140 


160 


6.2 ± 0.9 


6.4 ± 0.2 


123 ± 4 


D 


0.117 ± 0.028 


0.032 ± 0.006 


80 


100 


3.9 ± 0.8 


4.0 ± 0.2 


215 ± 12 


E 


0.100 ± 0.028 


0.042 ± 0.006 


80 


100 


3.8 ± 0.8 


4.4 ± 0.2 


148 ± 8 


F 


0.068 ± 0.025 


0.023 ± 0.016 


60 


80 


2.4 ± 0.8 


2.7 ± 0.2 


61 ± 5 


G 


0.064 ± 0.018 


0.007 ± 0.018 


40 


60 


1.9 ± 0.7 


1.4 ± 0.3 


106 ± 24 


H 


0.076 ± 0.019 


0.016 ± 0.020 


40 


60 


2.4 ± 0.8 


2.6 ± 0.3 


120 ± 15 


I 


0.052 ± 0.009 


0.022 ± 0.020 


25 


35 


1.9 ± 0.5 


1.1 ± 0.3 


174 ± 52 



